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ABSTRACT 

The  computer  storage  and  execution  time  requirements 
necessary  to  implement  a  detection  and  tracking  technique  based 
on  sequential  hypothesis  testing  have  been  examined.  Mathemati¬ 
cal  expressions  which  relate  these  requirements  to  sonar  system 
parameters  and  the  detection  algorithm  parameters  are  derived. 

It  is  shown  that  the  technique  can  be  implemented  on  modest, 
state-of-the-art  digital  computers.  In  addition,  future  studies 
designed  to  increase  the  shipboard  utility  of  the  detection  and 
tracking  technique  are  discussed. 
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1 .  INTRODUCTION 


1.1  PURPOSE 

The  formulation  of  a  computer  aided  detection  technique 
has  been  described  in  a  recent  technical  memorandum.^  This  tech¬ 
nique  is  based  on  sequential  hypothesis  testing.  The  technique 
employs  a  digital  computer  to  examine  in  real  time  the  output  of 
a  sonar  system,  reduce  the  volume  of  data,  perform  ping-to-ping 
Integration, and  generate  useful  displays.  A  question  of  immediate 
practical  interest  concerns  the  computer  storage  and  execution 
time  requirements  necessary  to  implement  this  technique  for  ship¬ 
board  usage.  The  purpose  of  this  memorandum  is  to  illustrate  how 
computer  requirements  can  be  estimated  and  how  the  sonar  system 
parameters  affect  the  computer  requirements. 

In  addition  to  these  practical  considerations,  a  dis¬ 
cussion  of  theoretical  problems  requiring  future  study  is  given. 
These  problem  areas  are  directly  related  to  shipboard  utilization 
of  the  likelihood  ratio  processor. 

1.2  BACKGROUND 

The  sequential  hypothesis  testing  procedure  considered 
in  this  note  is  discussed  in  Ref.  1.  Briefly  the  procedure  is  a 
method  of  deciding  between  two  hypotheses  Hq,  a  track  is  noise, 
and  Hp  a  track  is  signal-plus-noise.  The  quantity  tested  is  the 
logarithm  of  the  likelihood  ratio,  the  ratio  of  the  probability 
density  function  associated  with  to  the  probability  density 
function  associated  with  Hq.  Two  decision  thresholds  for  the 
logarithm  of  the  likelihood  ratio  are  derived  from  assigned 
probabilities  of  wrong  decisions:  a  lower  threshold  below  which 
it  is  assumed  that  the  track  is  noise  and  is  rejected,  and  an 
upper  threshold  above  which  it  is  assumed  that  a  track  is  a  tar¬ 
get  and  is  displayed.  For  samples  between  the  two  thresholds 
information  is  retained  and  combined  with  results  of  succeeding 
ping  cycles  until  the  joint  likelihood  ratio  crosses  the  upper 
or  lower  threshold. 
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The  data  flow,  shown  in  Fig.  I,  begins  with  the  outpxit 
of  the  sonar  processor  entering  a  Preliminary  Data  Rednction 
section  where,  after  thresholding,  a  single  ping  event  package 
containing  a  position  vector  and  the  log  likelihood  ratio  of  the 
event  is  created.  The  reduced  output  continues  to  the  New/Status 
Linkage  section  where  possible  linkages  with  previously  estab¬ 
lished  tracks  contained  in  the  statvis  file  are  considered.  For 
each  event  in  the  status  file  four  functional  quantities  are 
maintained:  a)  the  position  vector  in  the  proceeding  echo  cycle, 

b)  the  predicted  position  vector  in  the  present  echo  cycle,  c) 
a  position  variance  vector,  and  d)  the  log  likelihood  ratio  of 
the  event.  If  the  position  vector  of  a  single  ping  event  from 
the  Preliminary  Data  Reduction  section  lies  in  a  volvime  defined 
by  the  expected  position  vector  and  variance  vector  of  an  event 
in  the  status  file,  a  linkage  is  considered.  If  the  combined 
log  likelihood  ratio  is  greater  than  the  lower  decision  threshold, 
a  new  status  unit  is  created  and  stored  in  a  new  status  file. 

The  reduced  sonar  output  from  the  Preliminary  Data 
Reduction  also  goes  to  the  Secondary  Data  Reduction  section 
where  the  log  likelihood  ratio  of  the  single  ping  event  is 
tested.  If  this  is  above  the  lower  decision  threshold,  a  new 
status  unit  is  created  and  stored  in  a  new  status  file. 

The  status  file  from  the  previous  echo  cycle  is  pro¬ 
cessed  by  the  Status  Data  Reduction  section.  The  purpose  of  this 
reduction  is  to  maintain  well  established  tracks  even  thovjgh  on 
one  echo  the  new  data  does  not  pass  the  lower  threshold  test. 

The  log  likelihood  ratio  of  each  entry  in  the  status  file  is 
degraded  and,  if  it  is  still  above  the  lower  decision  threshold, 
a  new  entry  is  created  and  stored  in  a  new  status  file. 

The  final  process  is  carried  out  by  examining  the  log 
likelihood  ratios  of  the  newly  created  status  units.  Those  that 
exceed  the  upper  decision  threshold  arc  appropriately  displayed 
to  the  operator. 
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2 .  COMPUTEP  REQUIREMENTS 

2.1  INTRODUCTION 

Prime  considerations  in  implementing  a  likelihood  ratio 
processor  for  shipboard  usage  are  the  computer  storage  and  execu¬ 
tion  time  requirements.  These  requirements  are  determined  to  a 
large  extent  by  certain  design  parameters.  This  section  lists 
the  pertinent  parameters  and  relates  these  parameters  to  a  theo¬ 
retical  procedure  for  estimating  computer  requirements.  Numeri¬ 
cal  examples  are  included. 

2.2  DESIGN  PARAMETERS 

The  basic  design  parameters  which  should  be  taken  into 
consideration  when  implementing  the  computer  process  are: 

1.  The  output  statistics  of  the  sonar  signal  processor 
in  the  presence  of  noise  and  signal-plus-noise  and  the  output 
bandwidth  of  the  sonar  processor. 

2.  The  time  between  pings. 

3.  The  number  of  dimensions  with  which  possible  tar¬ 
gets  are  to  be  tracked  (range,  bearing,  Doppler,  etc.). 

4.  The  number  of  resolution  cells  in  each  dimension. 

5.  The  size  of  a  resolution  cell  in  each  dimension. 

6.  The  dynamic  constraints  on  the  rate  of  change  of 
each  dimension. 

7.  The  desired  display  clutter  rate. 

8.  The  minimum  signal-to-noise  ratio  of  interest. 

These  eight  items  will  be  divided  into  two  categories,  sonar 
parameters  and  decision  parameters,  and  are  discussed  in  the 
next  two  sections. 

2.2.1  Sonar  Parameters 

Since  the  sequential  testing  procedure  is  a  statisti¬ 
cal  test,  it  is  not  surprising  that  the  output  statistics  of  the 
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sonar  processor  are  the  most  important  factor  in  determining 
computer  requirements.  The  output  statistics  of  noise  alone  and 
signal-plus-noise  must  be  known  in  order  to  calculate  the  like¬ 
lihood  ratio  and  certain  probabilities  of  the  sequential  test 
continuing.  Although  it  is  easier  to  manipulate  the  probability 
density  functions  if  they  are  in  a  closed  form,  a  tabulated 
function  can  be  used  with  the  aid  of  a  digital  computer. 

The  output  bandwidth  of  the  sonar  is  a  factor  in 
determining,  among  other  things,  the  rate  of  independent  samples 
to  be  expected.  This,  in  turn,  has  a  bearing  on  the  number  of 
samples  to  be  allowed  for  each  range  or  time  resolution  cell. 

The  time  between  pings  when  compared  to  the  program 
execution  time  determines  whether  the  process  will  work  in  real 
time  or  not. 

2.2.2  Decision  and  Tracking  Parameters 

The  total  number  of  resolution  cells  is  of  prime  impor¬ 
tance.  This  number  is  given  by  the  product  of  the  number  of 
resolution  cells  allowed  in  each  dimension.  Some  possible  dimen¬ 
sions  are  range,  bearing  or  beam  number,  and  Doppler.  The  nvimber 
of  range  resolution  cells  is  given  by  the  time  period  searched, 
usually  the  time  between  pings,  divided  by  the  time  of  passage  of 
a  sonar  signal  through  a  range  resolution  cell. 

The  maximum  dynamic  range  allowed  in  any  dimension  is 
given  by  the  product  of  the  number  of  resolution  cells  in  that 
dimension  and  the  size  of  each  cell. 

Another  factor  which  affects  computer  loading  is  the 
maximum  rate  of  change  allowed  for  a  target's  movements  in  each 
dimension.  This  is  reflected  in  the  allowable  size  of  the 
variance  vector  in  each  dimension.  The  present  procedure  uses 
two  fixed  variances:  a  large  one  associated  with  single  ping 
events  allowing  for  any  reasonable  opening  or  closing  range  rate 
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and  a  smaller  one  associated  with  linked  tracks  allowing  for 
variations  in  the  observed  range  rate.  The  size  of  these  vari¬ 
ances  affects  the  number  of  linkages  to  be  made. 

The  most  important  result  obtained  by  this  process  is 
to  have  a  high  probability  of  detection  with  a  low  clutter  rate. 

The  clutter  rate  is  affected  by  the  specified  probabilities  of 
terminating  the  sequential  testing  procedure  with  an  incorrect 
decision  and  the  design  signal- to-noise  ratio  specified  to  deter¬ 
mine  the  log  likelihood  ratio. 

For  a  specified  minimum  detectable  signal-to-noise 
ratio,  a  design  signal-to-noise  ratio  should  be  chosen  such  that 
a  sample  of  the  minimum  signal-to-noise  ratio  will  be  assigned 
a  log  likelihood  ratio  value  of  zero.  If  this  is  done,  then 
there  will  be  a  .5  probability  of  detection  for  a  track  with 
that  average  signal-to-noise  ratio. 

2.3  AVERAGE  NUMBER  OF  STATUS  UNITS  REQUIRED 

The  objective  of  this  section  is  to  derive  a  formula 
for  the  average  number  of  status  units  created  in  one  ping  cycle 
by  the  sequential  testing  procedure  outlined  above.  It  is  assumed 
that  the  procedure  has  been  operating  in  a  noise  only  environment 
long  enough  to  reach  a  steady-state  condition,  that  is,  as  many 
status  units  are  being  created  as  discarded.  The  assumption  of 
noise  only  is  not  unreasonable  since  a  true  target  track  would 
produce  one  status  unit  and  possibly  a  few  status  units  describing 
branching  noise  tracks. 

By  considering  the  output  statistics  of  the  sonar  pro¬ 
cessor,  it  is  possible  to  find  the  probability,  Pj,  of  an  inde¬ 
pendent  sample  exceeding  a  specified  lower  threshold  T  corresponding 
to  the  one  in  the  Preliminary  Data  Reduction  section.  It  is  assumed 
that  the  maximum  independent  sample  in  a  resolution  cell  will  be 
chosen.  The  probability  that  at  least  one  of  the  Nj  independent 
samples  in  the  cell  will  exceed  the  threshold  T  is  given  by 
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Pi  *  1  -  (1  -  Pi)""'  .  (1) 

As  was  pointed  out  before,  the  total  number  of  resolu¬ 
tion  cells,  N^,  is  given  by  the  product  of  the  number  of  resolu¬ 
tion  cells  in  each  dimension. 

The  calculation  of  the  number  of  status  units  created 
or  continued  in  a  single  ping  will  be  divided  into  three  parts. 
The  first  source  is  from  the  Secondary  Data  Reduction  section. 
The  next  two  will  be  from  the  New/Status  Linkage  section,  but 
one  source  will  consider  only  the  linkage  of  single  ping ’events, 

and  the  other,  linkages  of  tracks  with  two  or  more  samples  or 
events . 

Another  possible  source  of  status  units  is  the  Status 
File  Reduction  section  of  the  processor.  This  section  is  mainly 
to  propagate  strong  tracks  that  have  perhaps  not  linked  in  the 
present  ping  cycle.  Since  only  noise  is  being  considered,  this 
section  will  not  contribute  many  status  units. 

For  any  sample  taken  from  a  resolution  cell,  a  like¬ 
lihood  ratio  L(x^)  is  formed  by 


where  p^(xp  is  the  probability  density  function  for  signal-plus- 
noise  and  Pq(x^)  is  the  probability  density  function  for  noise 
alone.  If  L(x^)  is  greater  than  B,  the  lower  decision  threshold 
a  status  unit  will  be  formed  by  the  Secondary  Data  Reduction 

section.  The  probability  that  the  status  unit  will  be  formed 
is  given  by 


P(B  <  L(x^))  =  J  Pq(x)  dx  , 


(3) 
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Pi  =  1 


1  -  P(B  <  L( 


xi)))  ' 


(A) 


The  probability  that  a  status  unit  will  not  be  formed 
is,  of  course,  1  -  Pi;  hence,  the  number  of  status  units  N’l 
formed  is  a  binominally  distributed  random  variable  with  mean 


^1  =  \ 


1  » 


(5) 


and  standard  deviation 


Ni  “  y^R  ^1  "  ^1^ 


(6) 


The  next  source  of  status  units  to  be  considered  is 
the  linkage  of  status  units  representing  single  ping  events  of 
the  previous  ping  cycle.  The  probability  P2  that  the  like¬ 
lihood  ratio  will  be  greater  than  the  lower  decision  threshold 
after  two  samples,  Xi,  X2,  have  been  considered  is  given  by 

P2  =  P[L(xi,X2)  >  B  I  Xi  >  T,  L(Xi)  >  B,  X2  >  T],  (7) 


or 

p  _  Is!  ‘^"'2  ‘^'^1 

2  “  p  p  » 

where  the  region  S  is  defined  by  the  conditions 
T  <  Xi^  <  • 

T  <  X2  <  ■ 

B  <  L(xi)  <  • 

B  <  L(xi,X2)  <  “  . 


(8) 
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On  the  average  the  number  N2  of  status  units  created 
by  single  ping  events  linking  with  noise  peaks  will  be 

N2  =  Ni  •  (Nl  P^)  •  .  (9) 

and  the  standard  deviation 


Nl  Pj 


P2(l  -  Pj) 


The  last  source  of  status  units  is  tracks  that  repre¬ 
sent  at  least  two  samples.  To  find  the  number  of  these  it  is 
necessary  to  know  how  many  samples,  on  the  average,  it  takes  to 
reach  a  decision  that  the  track  represents  noise.  The  direct 

computation  of  this  number  is  a  complex  procedure;  however,  by 

2 

certain  approximations  it  is  possible  to  find  a  good  value  for 
the  average  sample  number  of  all  possible  tests*.  Let  be 
the  average  number  of  samples  to  reach  a  decision  when  noise 
alone  is  present.  Formally  is  given  by 

m 

'I  1  •  Pi  ’  (ID 

i=l 

where  p^  is  the  probability  that  the  test  will  end  on  the  ith 
sample. 

If  N'  is  the  average  number  of  samples  for  tests  that 
have  more  than  two  samples,  then 


’"f^'he  probability  density  function  to  be  used  in  determining  the 
average  sample  number  is 

P'W  -  ^  PC*  <  X  I  T  <  X)  -  ^ 


p(x)/P(t  <  x) 


T  <  X 
otherwise 
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N 


F1P2 


av 


^  •  N;v  +  1  •  (1  -  p;^)  +  2  .  (pi)(l  -  P^),  (12) 


or 


^av  “  ^  2  • 


(13) 


By  simulating  the  sequential  testing  procedure  on  a  computer,  it 
is  possible  to  calculate  an  estimate  of  the  standard  deviation 
^av*  foe  the  examples  considered  below  the  estimate  was  on 
the  order  of  .5  when  was  in  the  range  of  3  to  4. 

The  average  number  of  status  units  which  have  been 
linked  in  established  tracks  can  be  expressed  as 

N3  -  N2  •  (Ng  .  P^)  .  (n;^  -  2)  ,  (14) 

where  Ng  is  the  number  of  resolution  cells  in  the  small  variance 
vector.  The  density  function  of  N'  ,  is  not  known:  hence,  the 
standard  deviation  of  N3  cannot  be  found.  It  seems  reasonable, 
however,  to  proceed  as  in  and  N2  and  assign  a  value  to  the 
standard  deviation  of  as  follows 


N. 


Pr) 


2 

av 


(15) 


2 

If  Nj^  and  are  large,  the  binomial  distributions 

may  be  approximated  by  Gaussian  distributions  with  the  given 
means  and  standard  deviations.  The  average  total  number  of 
status  units  for  steady-state  noise  only  conditions  is 


tot 


IT,  4-  IT.,  +  TTn 


(16) 


with  standard  deviation  of 
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C 


N 


tot 


(17) 


2.4  EXAMPLE 

The  example  to  be  considered  is  the  same  as  discussed 
in  Ref.  1  with  the  additional  consideration  of  a  multibeam  case. 
The  sonar  processor  consists  of  a  linear  correlator  followed  by 
an  envelope  detector.  The  bandwidth  is  100  Hz.  There  are  1500 
range  resolution  cells  and  in  the  multibeam  case,  five  beams. 

The  time  between  pings  is  assumed  to  be  15  sec.  The  large  fixed 
variance  for  single  ping  event  status  units  is  ±20  range  resolu¬ 
tion  cells  to  accommodate  targets  with  range  rates  of  ±20  knots 
and  smaller  range  variance  of  ±2  cells  to  allow  a  variation  of 
±2  knots  in  the  observed  range  rate.  In  the  multibeam  case  the 
beam  variance  is  ±1  in  both  cases.  In  determining  the  decision 
threshold  the  probability  of  accepting  a  noise  track  as  a  target 
is  assumed  to  be  0.001  and  that  of  rejecting  a  target  to  be  0.05. 
Since  only  normalized  Gaussian  noise  is  present,  the  output 
statistics  Pq(x)  are  Rayleigh: 

.x2 

where  N  is  the  rms  value  of  the  noise  before  envelope  detection. 
The  likelihood  ratio  is  given  in  Ref.  1  as 

L(x)  -  I^(x  S/N^)  exp(-S^/2N2)  .  (19) 

In  that  report  a  linear  approximation  of  the  logarithm  of  the 
likelihood  ratio  t(x)  is  derived  as 

•t(x)  -  ax  +  b,  (20) 
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where 

a  =  (S/N  -  l/6.74)/N  , 

b  -  3.87  S/N  -  ^(S/N)^  -  ^  log(7.75-  S/N)  -  a  •  3.87N  . 

Using  Eq.  (20)  it  is  relatively  easy  to  find  bounds  of  the  ran¬ 
dom  variables  when  compu*-ing  and  P2.  Since  each  resolution 
cell  will  contain  one  independent  sample. 


Pi  =  J  Po(^)  dx  , 


(21) 


logB-b 


and 


^2  = 


r 


r 


dx.  .  (22) 


logB-b  logB+logEpN,  -2b 

a  max( - - -  Xj^.T) 


In  the  case  considered  T  *  ;  hence  P_ 

a  ’  T 

beam  case  with  design  S/N  of  12  dB, 


Pj^.  For  the  single 


213 

120 


";v  '  3.5 


and  for  the  multibeam  case  , 
Pj^  -  .213 
P2  “  .0628 

"av  “  3.2 
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The  average  number  of  status  units  is  744  and 

7149,  respectively,  with  standard  deviations  of  23.3  and  58.4, 
respectively.  If  the  design  S/N  is  lowered  to  6  dB,  the  average 
number  becomes  2,481  and  20,070  with  standard  deviations  of 
35.6  and  72.3,  respectively. 

2.5  ACTUAL  STORAGE  REQUIREMENTS 

The  total  number  of  status  units  given  above  does  not 
give  the  true  computer  storage  requirements.  It  is  necessary 
to  find  the  number  of  computer  words  required  to  describe  each 
status  unit,  which  contains  a  position  vector  from  the  preceding 
echo  cycle,  the  expected  position  vector  in  the  next  echo  cycle, 
a  variance  vector,  and  the  joint  log  likelihood  ratios  of  the 
event . 

A  position  vector  can  include  range,  Doppler,  and 
bearing.  If  2048  or  less  range  resolution  cells  are  allowed, 
then  11  binary  bits  will  describe  the  range.  If  16  beam  and 
16  Doppler  cells  are  allowed,  then  4  bits  each  will  suffice. 

Hence  the  position  vector  from  the  preceding  echo  cycle  and  the 
expected  position  vector  for  the  next  echo  cycle  require  2(11+8)= 
38  bits.  The  variance  vector  associated  with  the  expected  posi¬ 
tion  vector  can  only  take  on  discrete  values;  therefore,  an 
indicator  of  3  bits  is  all  that  is  necessary  to  tell  the  program 
which  variance  vector  to  use.  The  likelihood  ratio  can  be  repre¬ 
sented  by  a  scaled  integer  number  of  12  bits,  giving  3  signifi¬ 
cant  decimal  digits. 

For  the  multibeam  case  with  Doppler  information  a 
total  of  53  bits  of  information  is  required  for  a  status  unit. 

If  the  fixed  word  length  of  a  computer,  like  the  AN/USQ-20A,  is 
30  bits,  the  first  word  could  contain  the  present  beam,  Doppler, 
and  range,  the  variance  indicator  and  the  predicted  beam,  and 
Doppler  all  presented  in  30  bits.  The  first  11  bits  of  the 
second  word  could  contain  the  predicted  range  and  the  last  12 


13 


6500  TRACOR  LANE.  AUSTIN.  TEXAS  78721 


bits  the  likelihood  ratio.  There  would  be  7  bits  unused;  hence 
a  total  of  two  computer  words  are  required  for  each  status  unit. 

The  process  as  it  is  now  programmed  requires  two  status 

file  storage  areas  of  equal  size,  one  for  the  previous  status 

file  and  one  for  the  file  being  created.  The  total  storage 

requirements  are  obtained  by  multiplying  the  total  number  of 

status  units  N.^.  by  4. 

tot 

The  numerical  results  of  the  example  given  show  that 
for  some  values  of  the  design  parameters  of  interest  the  storage 
requirements  can  be  met  by  a  modest,  state-of-the-art  computer. 

By  appropriate  changes  in  programming  and  storage  allocation  a 
wider  range  of  design  parameters  might  possibly  be  accommodated. 

2.6  AVERAGE  NUMBER  OF  COMPUTER  INSTRUCTIONS  REQUIRED 

If  the  storage  requirements  for  the  status  file  are 
within  reason  for  the  computer  being  considered,  the  next 
question  which  arises  is  the  expected  execution  time.  This 
depends  on  the  number  of  computer  instructions  to  be  carried 
out  during  a  ping  cycle  which  is,  in  turn,  related  to  the  num¬ 
ber  of  status  units  being  created.  To  find  the  number  of  instruc¬ 
tions  required,  the  processing  of  the  ping  cycle  will  be  separated 
into  three  parts,  the  Secondary  Data  Reduction,  the  New/Status 
Linkage,  and  the  Status  File  Reduction.  The  count  of  instruc¬ 
tions  will  vary  from  computer  to  computer  depending  on  the 
instructions  available.  The  numbers  quoted  here  are  based  on 
the  machine  language  listing  of  a  research  program  written  in 
UNIVAC  FORTRAN  V.  Speed  and  storage  were  sacrificed  where 
necessary  to  obtain  flexibility.  Careful  machine  language  pro¬ 
gramming  would  reduce  the  number  of  required  instructions. 

For  each  single  ping  event  out  of  the  Preliminary  Data 
Reduction  section  the  Secondary  Data  Reduction  section  requires 
instructions 
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II  =  li  +  (P^/Pt)Ii" 
=  9  +  (P,^/P^)22  , 


(23) 


where  Ij^'  is  the  nvimber  of  instructions  to  decide  whether  the 
single  ping  event  should  form  a  status  unit  or  not,  and  is 

the  number  of  instructions  to  create  a  status  unit. 


In  the  New/Status  Linkage  the  number  of  instructions 
required  for  each  sample  from  the  Preliminary  Data  Reduction 
section  is  approximately 


2Nt 

T  /  I  L  tot  -r  /  / 
^2  ^  ^2 


iTo  +  N, 

+  -i _ 1  1  '  " 

Vt  2 


(24) 


38  + 


2Nl  Ntot 


29  + 


N2  + 


N. 


80 


where  l2^  is  the  number  of  instructions  for  each  entry,  I2' ' 
the  number  required  to  consider  linkages,  ^-2' ' '  the  number 
required  to  create  a  new  status  unit. 

The  number  of  instructions  required  for  the  Status 
File  Reduction  section  is 


I 


3 


(25) 


where  I^'  is  the  number  of  instructions  to  consider  whether  to 
propagate  an  established  track.  Again,  the  numbers 
^2  *  ^2^  ’  ^2'  '  *  ^3^  based  on  the  requirements  of  the 

current  FORTRAN  program. 
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Since  the  total  number  of  input  samples  from  the  Pre¬ 
liminary  Data  Reduction  section  is  Ng  -  Nj^  P^,  the  total  number 
of  instructions  to  be  expected  for  each  ping  cycle  is 

^tot  "  ^T^^l  ^3  * 

If  Tj  is  the  execution  time  of  a  computer  instruction, 
the  total  execution  time,  is 


T 


tot 


tot 


(27) 


For  examples  given  above  using  a  Tj  of  8  microseconds,  the 
instruction  execution  time  for  an  AN/USQ-20A  computer,  the  number 
of  instructions  and  execution  times  for  various  design  param¬ 
eters  are 


Beam 

Design  S/N 

(dB) 

Ntot 

^tot 

Single 

12 

744 

408,284 

3.27 

Single 

6 

2,481 

2,331,636 

18.65 

Multi 

12 

7,149 

9,776,219 

78.21 

Multi 

6 

20,070 

50,736,557 

405.89 

The  above 

execution 

times  for  the  AN/USQ-20A 

do  not  look 

promising.  However,  it  should  be  emphasized  that  the  count  of  the 
number  of  instructions  to  execute  various  parts  of  the  processor 
is  based  on  a  machine  language  compilation  of  a  research-oriented 
FORTRAN  program.  When  adapting  the  process  to  a  particular  compu¬ 
ter  it  is  possible  to  streamline  the  coding  and  thereby  decrease 
the  program  execution  time.  Moreover,  modem  computers  have  greatly 
reduced  execution  times;  execution  times  of  one  microsecond  are  not 
uncomnon  with  present  generation  computers*.  If  the  execution  time 
for  one  Instruction  is  one  microsecond,  all  but  the  last  case  will 
run  in  real  time. 

*A  present  military  version  of  the  UNIVAC  1230  computer  has  a  core 
memory  storage  in  excess  of  30,000  words  and  an  execution  time  of 
1.8  microseconds. 
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3.  FURTHER  THEORETICAL  STUDIES 

3.1  LIMITING  COMPUTER  STORAGE  AND  EXECUTION  TIME 

The  previous  sections  show  that  lowering  the  design 
signal-to*noise  ratio  has  a  marked  tendency  to  Increase  the 
required  storage  and  execution  time,  and  It  Is  natural  to  ask 
how  these  requirements  might  be  limited  or  reduced.  One  or  a 
combination  of  the  following  study  areas  may  furnish  means  of 
reducing  computer  storage  and  execution  time  without  degrading 
the  detection  of  targets. 

1.  Adaptive  Adjustment  of  the  Threshold  in  the  Pre¬ 
liminary  Data  Reduction  Section.  This  would  allow  only  enough 
new  single  ping  events  into  the  sequential  processor  to  keep 
storage  near  the  maximum.  This  would  be  the  simplest  process 
to  implement,  but  it  introduces  the  possibility  that  small  tar¬ 
get  peaks  may  not  get  into  the  computer. 

2.  Periodically  Updating  the  Design  Signal-to-Noise 
Ratio.  An  increase  in  this  parameter  would  reduce  the  nximber  of 
status  units  in  the  status  file.  The  effect  of  this  on  the 
decision  process  must  be  studied  as  well  as  the  method  of  imple¬ 
menting  it. 

3.  Rejection  of  Multiple  Linkages.  The  present  pro¬ 
cessor  allows  more  than  one  track  to  share  the  same  peak  in  a 
ping  cycle,  each  track  having  a  different  predicted  position 
vector.  If  only  the  linkage  with  the  greatest  log  likelihood 
ratio  were  retained,  the  number  of  status  units  would  be 
decreased. 

4.  Linkages  Based  on  Other  Clues.  The  computer  pro¬ 
cessor  may  be  made  more  selective  by  applying  other  clues  to 
the  tracking  or  decision  algorithms.  For  example,  the  Doppler 
information  of  a  sample  could  be  compared  with  track  range  rate 
and,  if  consistent,  a  linkage  made.  If  aspect  information  is 
available,  then  it  could  be  applied. 
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5.  Modified  Log  Likelihood  Ratio.  An  interesting 
variation  to  the  sequential  test  would  be  to  choose  a  linear 
function  of  the  sonar  output  that  would  give  a  specified  average 
number  of  samples  for  a  decision  in  the  presence  of  noise  and 
for  a  signal  with  the  design  signal-to-noise  ratio.  Since  the 
likelihood  ratio  sequential  test  is  optimum,  any  variation  from 
it  would  be  sub-optimum  in  the  overall  system;  however,  the 
degradation  may  be  small  and  the  benefits,  faster  rejection  of 
noise  and  integration  of  signal-plus-noise  target  tracks,  could 
outweigh  the  degradation. 

3  4 

6.  Adaptive  Testing.  Two  methods  ’  of  adaptive 
sequential  testing  have  been  developed  recently.  While  not 
limiting  computer  storage,  they  would  allow  the  computer  pro¬ 
cessor  to  adapt  to  the  changing  conditions  of  the  target  or  the 
sea.  The  feasibility  of  implementing  these  methods  should 
receive  attention. 

3.2  SHIPBOARD  CONSIDERATIONS 

The  above  areas  explore  ways  of  improving  the  internal 
processes  of  the  sequential  detector.  There  are  other  items 
which  must  be  considered  before  a  useful  system  can  be  imple¬ 
mented  on  a  ship.  These  are: 

1.  Tracking  Procedures.  A  realistic  system  must 
adjust  the  predicted  position  vectors  contained  in  the  status 
units  when  the  ship  is  changing  course  and  speed.  The  manner 
in  which  this  is  done  depends  on  the  coordinate  system  used  by 
the  sequential  detection  process. 

2.  Additional  Target  Information.  The  computer 
should  certainly  be  able  to  handle  duties  other  than  detection. 
The  position  vectors  of  a  status  unit  for  the  last  several  pings 
should  be  available  for  display  if  the  log  likelihood  ratio  is 
large  enough.  Also,  range  rate,  Doppler,  course,  and  speed 
should  be  available  for  display  purposes.  Storage  of  multiple 
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ping  information  will,  of  course,  increase  computer  require¬ 
ments  . 

3.  Clutter  Rate.  The  display  clutter  rate  and  the 
probability  of  detection  as  a  function  of  the  design  parameters 
should  be  studied.  The  display  threshold  will  be  one  parameter 
studied. 

3.3  AREAS  TO  BE  STUDIED  IN  THE  IMMEDIATE  FUTURE 

It  is  doubtful  that  all  the  areas  outlined  above  can 
be  accomplished  under  the  present  contract.  Emphasis  will  be 
placed  on  the  following  areas. 

1.  Adaptive  Adjustment  of  the  Threshold  in  the  Pre¬ 
liminary  Data  Reduction.  The  threshold  adjustment  will  be  based 
on  the  number  of  status  units  created  in  a  ping  cycle.  The 
threshold  may  be  adjusted  once  each  cycle  or  periodically  during 
the  cycle.  The  effect  of  this  adjustment  on  detection  and  the 
algorithm  for  the  adjustment  will  be  studied. 

2.  Tracking  Procedures.  The  most  promising  tracking 
coordinate  system  appears  to  be  polar  coordinates  relative  to 
own  ship's  course.  This  system  can  be  modified  to  take  into 
account  own  ship's  maneuvers  by  adding  a  correction  vector  to 
each  position  vector  in  the  status  file.  How  often  this  updating 
is  necessary  and  exactly  how  to  calculate  the  correction  vector 
are  subjects  for  further  study. 

3.  Linkages  Based  on  Doppler.  The  most  promising 
way  to  make  the  computer  processor  more  selective  in  retaining 
linkages  would  be  to  use  Doppler  information  to  decide 

the  observed  range  rate  is  consistent  with  the  Doppler  informa¬ 
tion.  Doppler  information  would  also  allow  range  rate  informa¬ 
tion  to  be  applied  to  single  ping  events,  thus  the  window  widths 
(variances)  may  be  decreased  in  this  case.  Both  the  probability 
density  functions  involved  and  the  implementation  algorithm  will 
be  studied  in  detail. 
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4 .  PERFORMANCE  EVALUATION 


Once  the  major  practical  difficulties  relating  to  the 
feasibility  of  shipboard  implementation  have  been  resolved,  it 
would  be  most  desirable  to  conduct  a  thorough  performance  eval¬ 
uation  of  the  sequential  likelihood  ratio  processor.  It  is 
generally  felt  that  sequential  likelihood  ratio  processing  offers 
certain  advantages  over  standard  signal  processing  techniques 
such  as  single-ping  matched  filtering.  Two  of  these  advantages 
are  a  reduction  in  the  amount  of  extraneous  data  presented  to 
the  operator,  and  ping-to-ping  integration.  In  order  to  demon¬ 
strate  that  these  advantages  translate  into  improved  system  per¬ 
formance,  the  following  tests  should  be  conducted. 

The  performance  of  a  typical  active  sonar  detection 
system  should  be  measured  in  terms  of  receiver  operating  char¬ 
acteristic  (ROC)  curves.  To  do  this,  typical  ping  cycles  would 
be  synthesized  and  applied  to  a  signal  processor  such  as  a 
linear  correlator.  The  correlator  output  would  then  be  thresh- 
olded  and  applied  to  a  display  such  as  the  simulated  AN/SQS-26 
B-scan.  (This  display  has  been  simulated  in  the  TRACOR  display 
facility  for  an  AN/SQS-26  mutual  interference  study.)  By  con¬ 
ducting  observer  tests  with  the  display  it  would  be  possible  to 
obtain  ROC  curves  for  the  simulated  detection  system. 

Next,  a  sequential  likelihood  ratio  processor  imple¬ 
mented  on  a  shipboard  type  computer  such  as  the  AN/USQ-20  (or  a 
simulated  version  thereof)  should  be  inserted  between  the  signal 
processor  and  the  display.  Observer  tests  would  be  conducted 
and  new  ROC  curves  obtained.  In  this  way,  the  gains  in  overall 
system  performance  arising  from  the  use  of  sequential  likelihood 
ratio  processing  could  be  measured  in  a  quantitative  fashion. 
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5.  CONCLUSIONS 

Expressions  which  relate  computer  storage  requirements 
and  program  execution  time  to  parameters  of  a  sonar  system  and 
the  likelihood  ratio  processor  have  been  derived.  The  numerical 
results  indicate  the  processing  technique  can  be  implemented  on 
a  modest,  state-of-the-art  digital  computer.  For  less  advanced 
computers  such  as  the  AN/USQ-20,  the  likelihood  ratio  processor 
can  be  implemented  with  some  restrictions  imposed  on  the  pro¬ 
cessor  decision  and  tracking  parameters.  Theoretical  analyses 
for  the  immediate  future  will  be  aimed  at  modifications  to  the 
existing  algorithms  such  that  computer  requirements  will  be 
minimized. 
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